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of O I, N I, and N ll were derived from spectra obtained by the Far Ultraviolet 
Spectroscopic Explorer (FUSE), the column density of Dl was derived from a 
combination of our FUSE spectra and an archival HST GHRS spectrum, and 
the column density of Hi was derived from a combination of the GHRS spectrum 
and values derived from EUVE data obtained from the literature. We find the 
following abundance ratios (with 2cr uncertainties): Dl/Hl = (1.66±0.28) x 10~^, 
Oi/Hi = (3.63 ± 0.84) x 10"^, and Ni/Hi = (3.80 ± 0.74) x 10"^ The Nil 
column density was slightly greater than that of N I, indicating that ionization 
corrections are important when deriving nitrogen abundances. Other interstellar 
species detected along the fine of sight were C ll, C ill, O VI, Si ll, Ar I, Mg ll, and 
Fell; an upper limit was determined for Nlll. No elements other than Hi were 
detected in the stellar photosphere. 

Subject headings: ISM: abundances — ISM: clouds — stars: individual (HZ 43A) — 
stars: white dwarfs 



1. Introduction 

Deuterium is one of the primary products of big bang nucleosynthesis (BBN), and its 
primordial abundance provides a sensitive measure of the baryon density of the universe 
(Schramm & Turner 1998; Buries 2000). However, deuterium is consumed in stars far more 
quickly than it is produced by the P-P cycle or by any other known stellar nucleosynthetic 
process (Epstein et al. 1976); thus, measurements of the deuterium abundance at the present 
epoch can only provide lower limits to the primordial abundance. Measurements obtained 
at a variety of redshifts will sample the chemical evolution of the universe, and may permit 
extrapolation to a value for the primordial abundance of deuterium. Measurements of the 
abundance of deuterium and various products of stellar nucleosynthesis, such as oxygen, 
within our Galaxy will improve our understanding of chemical evolution as well as provide a 
lower limit to the primordial deuterium abundance. Recent reviews of deuterium abundance 
measurements can be found in Lemoine et al. (1999), Linsky (1998), and Moos et al. (2001). 

Linsky (1998) reviews measurements of the deuterium abundance for material in the 
local interstellar cloud (LlC), which are consistent with a common value for D/H of (1.5 ± 
0.1) X 10^^. The LlC and other nearby clouds are thought to be embedded in a large bubble of 
hot, low density gas (T IQi^K, n ^ 0.005 cm~^) that was probably produced by supernovae 
and stellar winds arising in the Scorpius-Centaurus OB association (Cox & Reynolds 1987; 
Frisch 1995). If this picture is correct, then the local interstellar medium (LISM) may contain 
an inhomogeneous mixture of clumps of older material swept up by the expanding bubble 
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along with the more-recently processed material in the bubble itself. Thus, even in this local 
environment one may find regions of gas with somewhat different evolutionary histories and 
different compositions. The same processes have been at work throughout the history of the 
Galaxy; hence measurements of deuterium and other abundances along numerous lines of 
sight will help determine not only the chemical evolution of the Galaxy, but also the relative 
timescales for chemical evolution and the mixing of material. 

In this paper we present an analysis of the line of sight to the nearby white dwarf 
HZ43A. Companion papers will present similar analyses of the lines of sight to G 191-B2B 
(Lemoine et al. 2001), WD 0621-376 (Lchncr ct al. 2001), WD 1634-573 (Wood ct al. 2001), 
WD 2211-495 (Hcbrard et al. 2001), BD +28° 4211 (Sonneborn et al. 2001), and FeigeUO 
(Friedman et al. 2001). An overview will be provided by Moos et al. (2001). 

High-resolution spectra of the DA white dwarf HZ43A, covering the far ultraviolet 
(FUV) wavelength range 905-1187A were obtained with the Far Ultraviolet Spectroscopic 
Explorer (FUSE) for the purpose of studying the deuterium abundance of the local interstellar 
medium (ISM). This fine of sight is promising for a study of D/H for several reasons: the star 
itself is moderately bright, exhibits a nearly featureless continuum, and has been accurately 
modelled; the absorbing interstellar gas appears to have only a single velocity component 
and the column density is very low, so many absorption lines that arc ordinarily on the flat 
part of the curve of growth are either unsaturated or exhibit only mild saturation effects. In 
addition, the H I column density along the line of sight can be determined by two independent 
methods: from the shape of the Lyman continuum observed by the Extreme Ultraviolet 
Explorer (EUVE), and from the Lyman-ct profile observed with the Goddard High Resolution 
Spectrograph (GHRS) onboard the Hubble Space Telescope {HST). Systematic uncertainties 
in determination of Hi column densities are often the limiting factor in measurements of 
the D/H ratio, so the availability of multiple independent measurements of N(H) is quite 
valuable. 

In Section 2 we discuss the line of sight to HZ43A and the properties of the star; 
in Section 3 we describe the observations and data reduction procedures; in Section 4 we 
describe our analysis procedures, our methods for minimizing systematic errors, and our 
measured column densities; and in Section 5 we summarize our results. 

2. Line of Sight and Stellar Properties 

The hot DA white dwarf HZ43A (WD 1314-^293) is a member of a non-interacting 
binary system with a dMe star. The cool secondary star has no significant UV fiux. The 
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coordinates of the star and some other basic characteristics are given in Table 1. The Galactic 
coordinates, / = 54.10°, b = +84.16°, are close to the north Galactic pole and the line of sight 
exits the LIC after traversing a very short distance: ~0.1pc (Redfield & Linsky 2000). The 
mean Hi density in the LIC is 0.1 cm~^ (Linsky et al. 2000), so the corresponding column 
density is only 3.1 x lO^^cm"^. The velocity of the LIC projected onto the HZ43A line 
of sight is -8.9km s~^ (all velocities given in this paper are heliocentric). The line of sight 
may intersect the "G" cloud that hes between us and the Galactic center, and the "NGP" 
cloud that hes between us and the North Galactic Pole (see Sfeir et al. 1999 and the review 
by Prisch 1995 for a discussion of the geometry of the nearby clouds lying within the local 
bubble). The boundaries of the clouds are not known in detail, but the identity of each 
absorber can be determined by its relative velocity. In particular, the velocity of the G cloud 
projected onto the HZ43A line of sight is -12.1km s^^. The distance to HZ43A has been 
determined by several methods to be approximately 65pc. Vennes et al. (1997) obtained 
67pc from a comparison of its apparent and absolute V magnitude (derived from a model 
atmosphere), while Dupuis et al. (1998) obtained a somewhat smaller value, 56-6 Ipc using 
a similar method. Ground-based parallax measurements have given 48-91pc (Margon et al. 
1976), 53-83pc (Dahn et al. 1982), and 55-81pc (van Altena et al. 1995). The only discrepant 
measurement, 25 - 44pc, is from the Hipparcos catalog (Clark & Dolan 1999); the origin of 
this discrepancy is not understood. The apparent line-of-sight velocity of the star (including 
gravitational redshift) was measured by (Reid 1996) to be -1-20.6 km s~^. 

The atmospheric parameters of HZ 43A are difficult to measure at visible wavelengths, 
because of the presence of the bright M star companion only ~ 3" away. Napiwotzki et al. 

(1993, N93 hereafter) and Finley et al. (1997, FKB hereafter) obtained optical spectra of 
both components and corrected the observed flux of the white dwarf by subtracting the 
contribution from the M star. N93 obtained T^s = 49,000 K and log (7 = 7.7 using NLTE 
models, while FKB determined T^s = 50,800 K and logg = 7.99 using LTE models. Ac- 
cording to Napiwotzki et al. (1999), the difference between both studies may illustrate the 
correction in T^s and log 51 that one must apply to transform the atmospheric parameters 
determined by LTE models to NLTE models. However, if we apply the Napiwotzki et al. 
(1999) correction (Aloggi = —0.13) to PKB's value (7.99), we obtain a larger gravity than 
N93 (7.7). 

Other spectral bands have been used to evaluate the effective temperature and gravity 
of HZ43A. For example, Holberg et al. (1986) fitted the lUELjman-a line profile, obtaining 
a larger temperature and gravity than the optical observations performed by N93 and FKB. 
Both Vennes (1992) and FKB indicated that the visual magnitude estimated by Holberg 
et al. (1986) was too faint. The visual magnitude is quite difficult to determine from the 
ground because of contamination by the companion star. If Holberg et al. (1986) had used 
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the magnitude determined from HST data by Bohlin et al. (1995), the effective temperature 
would have been Tefj = 50, 850 K, very close to that determined later by FKB. The complete 
Lyman series of HZ43A was analyzed by Dupuis et al. (1998) using ORFEUS and by Kruk 
et al. (1997) with the Astro-1 Hopkins Ultraviolet Telescope, while the EUV wavelength 
range was studied by Dupuis et al. (1995) and Barstow et al. (1995) using EUVE. If we take 
the average of the effective temperature and gravity for all measurements made from EUV, 
FUV, and optical observations, we obtain T^q — 50,300 K and logg^ = 7.88, with standard 
deviations of 700 K and 0.15 dex (see, e.g.. Table 1 of Dupuis et al. 1998). 

For the purpose of this study, we adopted the following atmospheric parameters: Tefj = 
50,900 K, and \ogg = 8.0. These are the pre-publication values of FKB used by Kruk 
et al. (1999) in the final Astro-2 calibration of the Hopkins Ultraviolet Telescope, and used 
presently in the FC/5'£^ calibration (the present FUSE Hux cahbration is defined by HZ43A, 
G191-B2B, and GD246). 

HZ 43A is unusual in that its photosphere appears to consist entirely of hydrogen, though 
stars with similar atmospheric parameters show traces of heavy elements. Barstow et al. 
(1995) set an upper limit on the helium abundance of 10~^, based on the EUVE spectrum 
of this star, and Dupuis et al. (1998) used an ORFEUS spectrum to set abundance limits 
of 10"^ for S, 10-s for C and N, and 10-^-^ for Si, P, and CI. Venues et al. (1991), N93, 
and Dupuis et al. (1995) came to the same conclusion, though with somewhat higher upper 
limits. Therefore, photospheric models for this star employing a pure hydrogen composition 
should give accurate results. 



3. Observations 

3.1. FUSE Observations 

FUSE comprises four independent co-aligned telescopes and spectrographs. Two of the 
channels contain optics coated with aluminum and overcoated with LiF, and the other two 
channels contain optics coated with SiC. Together the four channels span the wavelength 
range 905-1 187A. Most wavelengths in this band are sampled by at least 2 channels, and 
roughly one-third of the bandpass is sampled by 3 or 4 channels (lOOOA < A < IO8OA). 
Comparison of features in the spectra obtained in different channels can thus be used to test 
for the effects of fixed-pattern noise, backgrounds, or other instrumental characteristics that 
vary from one channel to another. Each spectrograph has 3 entrance apertures: LWRS (30" x 
30"), MDRS (4"x 20"), and HIRS (1.25" x 20"). Each aperture illuminates different regions 
of the detector and produces spectra with slightly different line-spread functions (LSFs). 
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Thus, if a source is observed through multiple apertures one can perform additional tests 
for the presence of fixed-pattern noise introduced by the detector, and one can determine 
if contamination by geocoronal emission is significant. Further information on the FUSE 
satellite and instrument on-orbit characteristics are provided by Moos et al. (2000) and 
Sahnow et al. (2000). 

Three FUSE observations of HZ 43A have been obtained, two in the LWRS aperture 
and one in the MDRS aperture. Basic information for the datascts is given in Table 2. 
The data were obtained in "histogram" mode, in which a two-dimensional spectral image is 
accumulated and downlinked. The LWRS observations will be described here for the sake of 
completeness: the signal-to-noise of the LWRS observations was too low in comparison with 
the MDRS data for them to significantly constrain the results. They were useful primarily 
for tests of fixed-pattern noise effects. Exposure durations for the LWRS observations ranged 
from just under 500 s to just over 800 s. The exposures in the MDRS observation were all 
approximately 485 s, apart from a few exposures that were shortened to roughly 200 s by de- 
layed target acquisitions. The MDRS observation was split into 4 intervals of approximately 
equal exposure time, each with a different X-position of the Focal Plane Assemblies (FPAs). 
The shift in FPA position causes a corresponding shift in the position of the spectrum on 
the detector along the dispersion direction. The 4 positions were at X = -57 /xm, 0/xm, 
+33 iim, and -M25 /xm with respect to the nominal position, corresponding to spectral shifts 
of roughly -12, 0, -|-5, and -|-19 pixels. These shifts span a range several times larger than the 
width of typical interstellar absorption lines, and several times larger than the width of the 
high-frequency components of the detector fixed-pattern noise. This so-called "focal-plane 
split" ("FP-split") procedure results in a considerable improvement in the signal-to-noise 
ratio achieved, even if the the spectra from the different FPA positions are simply shifted 
and added with no attempt to explicitly solve for the detector fixed-pattern noise. 

The data were processed using the most recent version of the FUSE calibration pipeline, 
CALFUSE version 1.8.7. Because there is no information available for the arrival time of 
individual photons in histogram mode, the pipeline applies exposure-averaged corrections for 
the Doppler shift induced by the spacecraft orbital motion and for the rotation of the gratings 
on orbital timescales. The primary mirrors also rotate slightly on orbital timescales, which 
can lead to significant changes in the position of the image of the source at the spectrograph 
entrance aperture. As a result, the effective spectral resolution can vary from one exposure 
to another and the zero-point of the wavelength scale varies from one exposure to another. 

The rotations of the primary mirrors during an exposure can cause the stellar image 
to fall outside of the corresponding spectrograph aperture, resulting in loss of signal in that 
channel for all or part of an exposure. The Fine Error Sensor presently used to provide 
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pointing information to the satellite attitude control system is part of the optical path 
in the LiFl channel; hence the target is always maintained at a fixed position in the LiFl 
spectrograph aperture. The alignment of the other channels is usually maintained so that the 
rotations of the primary mirrors occurring on orbital timescales do not shift the source image 
at each focal plane out of the corresponding LWRS aperture. For the M1010501 observation, 
no flux was lost in any of exposures; for P1042301 no flux was lost except that the SiCl 
channel fluxes were low by 7%, 10%, and 19% in exposures 13, 16, and 20, respectively. The 
SiCl spectra obtained in these exposures were scaled accordingly before being combined with 
the other exposures. For MDRS observations, peakups are performed twice in each orbit 
in order to try to correct for the mirror rotations by moving the spectrograph apertures 
(FPA)s in the dispersion direction to maximize the signal prior to the start of an exposure. 
For the PI 042302 MDRS observation, this procedure was fairly successful, providing effective 
exposure times of 82% for LiF2, 62% for SiCl, and 76% for SiC2. The necessity of moving 
the FPAs to compensate for rotations of the primary mirrors means that there is actually 
a distribution of FPA positions clustered about the mean offsets that define the FP-split 
sequence described above. 

The spectra obtained in different exposures generally need to be co-aligned (shifted to 
a common wavelength zero-point) prior to combining. This is caused in part by the lack of 
corrections for primary mirror rotations in the present pipeline, and in part by shortcomings 
in the present corrections for grating rotations. For HZ 43A there are no narrow photospheric 
lines and very few strong interstellar absorption lines. For the LiF channels, we fit a Gaussian 
profile to the interstellar C ll 1036. 337A line in each exposure and shifted each exposure to 
locate this line at a common velocity. For the SiC channels, we performed a cross-correlation 
on the the interstellar Lyman lines between 915A and 930A (fitting the ClI 1036. 337A line 
gave essentially identical results). In order to have enough signal to measure line positions, 
we discarded spectra for which the fiux was below half of the nominal value for its channel. 



3.2. GHRS Observations 

Archival HS T /GRRS spectra were used to determine the velocity structure of the in- 
tervening gas and to measure Lyman-a absorption of both Hi and D I. A log of the GHRS 
observations is given in Table 3. All of the data were acquired after the installation of the 
Corrective Optics Space Telescope Axial Replacement (COSTAR) unit. The cchcllc-mode 
data, all of which were acquired through the large science aperture (LSA; 1.74"x 1.74"), 
have a fine spread function with a narrow core of resolution 3.5 km s~^ (FWHM) and broad 
wings (~15 km FWHM) with relatively low power (see Appendix B of Howk et al. 1999). 
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Our reduction of these data follows that of Howk et al. (1999). The basic calibration 
makes use of the standard CALHRS pipeline with the best calibration reference files avail- 
able as of the end of the GHRS mission. The CALHRS processing includes conversion of 
raw counts to count rates and corrections for particle radiation contamination, dark counts, 
known diode nonuniformities, paired pulse events and scattered light. The wavelength cal- 
ibration was derived from the standard calibration tables based upon the recorded grating 
carrousel position. The absolute wavelength scale should be accurate to ~1 resolution ele- 
ment, or ~3.5 km s"^ for the echelle-mode data used in this work. 

All of the cchcllc-modc observations employed here were taken at four separate grating 
carrousel positions to place the spectrum at different positions along the 500-diode science 
array of the GHRS. This is done to mitigate the effects of the fixed-pattern noise that can 
limit the signal-to-noise ratio and produce artifacts in GHRS spectra. We have co-added each 
of the individual sub-exposures by aligning the interstellar hnes of interest (e.g., Dl, Mgll, 
and Fell) in each of the observations. We used the algorithm of Cardelli & Ebbets (1994) 
to explicitly solve for the fixed-pattern noise spectrum in each of the datasets. However, we 
found it unnecessary to apply this correction because the signal-to-noise ratio of the data did 
not necessitate the removal of the low-power features and no prominent features were found 
(see discussion in Howk et al. 2000). We note that the coaddition of the shifted datasets 
reduces the strength of any fixed-pattern features by a factor of four. 

The inter-order scattered light removal in GHRS cchcllc-modc data discussed by Cardelli 
et al. (1990, 1993) is based upon extensive prc-flight and in-orbit analysis of GHRS data and 
is used by the CALHRS routine. The coefficients derived by these authors are appropriate for 
observations made through the small science aperture (SSA). Corrections to these coefficients 
are often required to bring the cores of strongly saturated hnes observed through the LSA 
to the correct zero level. Typically this is done by adjusting the spectrum zero level by 
an amount equal to a fraction of the average net flux in each order. This "d" coefficient 
(Cardelli et al. 1990, 1993) accounts for the residual effects of the local echelle scattering. 

For the Ech-B data covering Mg ll and Fe ll, no strongly saturated hnes were available for 
evaluating the quality of the background removal, so we adopted the standard d-coefficients 
from the CALHRS processing. The Lyman-a absorption profile for HZ 43A should be flat- 
bottomed at zero flux near line center. Although the presence of strong Lyman-a airglow 
confuses the matter somewhat, it is clear that the use of the standard d-cocfficient (d=0.052) 
over- subtracts the core of the Lyman-a profile. We have used the observed flux level over 
the wavelength range 1215.74 to 1215. 83A, assuming this represents the true zero point, to 
adjust the d-coefficient. The adopted d-coefficient is d = -0.006±0.033. 

The G160M observations of the 1192-1228A range were acquired through the Small Sci- 
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ence Aperture (0.22" x 0.22"; SSA). These data have a resolution of ~19.6km s"^ (FWHM). 
The scattered hght properties of this holographicahy-ruled first-order grating are exceUent. 
No corrections to the CALHRS background subtractions were necessary. 

4. Analysis 
4.1. Method 

We attempted to determine column densities by more than one method whenever possi- 
ble. The two primary methods were profile fitting and curve-of-growth analyses. These two 
methods differ in their sensitivity to some of the sytematic uncertainties that affect column 
density measurements. A curve-of-growth analysis was not feasible for Hi, as essentially all 
of the lines other than Lyman-a are on the flat part of the curve of growth. Instead, Hi 
column density measurements based on fits to the EUV spectrum of HZ 43A were taken from 
the literature to provide a consistency check on our results. 

Profile fitting was used to determine the column densities of Hi, Dl, Ol, and Nl. The 
profile fitting was performed using the program Owens. f, written by one of us (M.L.). This 
program uses minimization to adjust temperatures, turbulent velocities, and relative 
velocities for multiple gas components, to adjust column densities for each species present in 
each gas component, and to adjust polynomial fits to continuum fluxes, line spread functions, 
background levels, etc. for each specified spectral region. Each quantity can be adjusted by 
the program or held fixed to user-specified values. Further information on this program is 
provided in the companion paper by Hebrard et al. (2001). 

The FUSE and GHRS spectra were normalized prior to fitting by means of an NLTE 
stellar model. The continuum in the vicinity of individual absorption lines was then fit 
by a low-order polynomial as required to account for any residual flux calibration errors 
or flat-field effects. We use the program TLUSTY (Hubeny & Lanz 1995) to calculate the 
atmospheric structure of HZ 43A under the assumption of NLTE, and then we compute the 
detailed emergent fiux using the program SYNSPEC (Hubeny 2000, private communication). 
The hydrogen line profiles are obtained from the Stark broadening tables computed by Lemke 
(1997) using the program of Vidal et al. (1973). We use a NLTE model atmosphere to 
calculate the emergent flux, because NLTE effects are important in the core of the Lya 
line, which is formed high in the atmosphere where departure from LTE is significant (see, 
e.g., Wesemael et al. 1980). The model parameters used were as discussed in Section 2: the 
composition was pure hydrogen, Tcff= 50900K, and logg= 8.00. The model spectrum was 
then normalized to V = 12.914 (Bohlin et al. 1995). Except for NLTE effects in the narrow 
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cores of the Lyman lines it is essentially identical to the LTE model described in detail by 
Kruk et al. (1999) and references therein. 

The Lyman-a profile was also fit using a different program written by one of us (J.C.H), 
and the profile fitting code of Ed Fitzpatrick (Spitzer & Fitzpatrick 1993) was used to provide 
a consistency check for certain test cases. 

We have derived column densities for N I and O I by fitting the measured equivalent 
widths with a single-component Gaussian (Maxwellian) curve of growth (Spitzer 1978; Jenk- 
ins 1986). The analysis was performed using software adapted from that described by Savage 
et al. (1990), which varies the Doppler parameter, b, and column density so as to minimize 
between the observed equivalent widths and a model curve of growth. The stellar continuum 
in the vicinity of each line was estimated using a low-order Legendre polynomial fit to the 
data. The dispersion of the data about the fit was used to estimate the uncertainties in the 
final integrated equivalent width measurements as discussed by Sembach & Savage (1992). 
This estimate includes contributions from both Poisson noise and uncorrected high-frequency 
fixed-pattern detector features. The measurement uncertainties quoted include contributions 
from this noise estimate, uncertainties in the Legendre fit parameters (Sembach & Savage 
1992), and estimates of the systematic uncertainties inherent in the continuum placement 
and velocity integration range. Atomic data for both the profile fitting and the curve of 
growth analysis were taken from Morton (1991, Morton 2001, private communication). 

The integrated equivalent widths for several interstellar species are presented in Table 
4, and the results from the curve-of-growth fits are given in Table 5. For most species other 
than N I and O I the number of available transitions was insufficient to fit a curve of growth 
and determine a b-value; instead, the measured equivalent width was combined with a curve 
of growth computed using the O I b-valuc to derive the column density. For these cases 
the table also lists the difference between the column density derived in this manner and 
that derived using a straight integration of the apparent column density profiles (Savage & 
Sembach 1991), providing a measure of the saturation correction. 

4.2. Systematic Uncertainties 

Column density determinations are subject to numerous possible systematic uncertain- 
ties. These are described in the following section, along with brief descriptions of how these 
were addressed. The magnitude of each effect on the column densities will then be addressed 
in the discussion of each species as appropriate. 
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4-2.1. Overlapping velocity components along the line of sight 

One of the major difficulties typically encountered in measuring column densities is 
accounting for overlapping absorption profiles resulting from clouds with similar velocities. 
As we will show below, there appears to be only a single absorption component along the 
line of sight to HZ43A, so this complicating factor is not present. 

4.2.2. Narrow absorption features in the stellar spectrum 

Stellar photospheres ordinarily contain trace quantities of numerous elements that can 
give rise to narrow absorption features in the emergent spectrum. In some cases the stellar 
spectrum is well-understood and the presence of such narrow features can be predicted 
and possibly even modelled. In other cases, the wavelengths of the various transitions are 
poorly known and it is impossible to be certain how much of an absorption profile is due to 
interstellar gas and how much results from the photospheric spectrum. However, as discussed 
in Section 2, the atmosphere of HZ 43A appears to be composed solely of hydrogen. Thus, 
this potential complicating factor is also not a concern for our analysis of HZ43A. 

4.2.3. Uncertainties in the stellar Lyman-a profile 

The high gravity of HZ43A causes the photospheric Lyman fines to be quite broad, 
hence the stellar spectrum is varying slowly in the vicinity of all of the interstellar lines being 
fit other than Hi Lyman-a. Because we are fitting the damping wings of the interstellar 
Hi Lyman-a profile, spanning roughly 6A, we need to be able to model the photospheric 
absorption very accurately. In addition, the Hi column density along the line of sight 
to HZ43A is sufficiently low, and the relative velocity of the star and interstellar gas is 
sufficiently high, that the NLTE core of the photospheric Lyman-a profile may affect the 
red wing of the combined profile. HZ 43A has been studied in great detail, however, and 
model atmospheres and synthetic spectra of HZ43A, along with those of several similarly 
well-studied DA white dwarfs, are used to define the fiux calibration for most orbiting UV 
spectrometers (Bohfin 2000, and references therein). The parameters describing HZ43A are 
discussed in Section 2, and the fits to the Lyman-a profile are discussed in Section 4.3.2. 



4 .2. 4- Continuum placement errors 



Errors in placement of the continuum or the background will clearly affect the column 
density determination. In most regions of the spectrum the continuum appears to be smooth 
and slowly-varying in the vicinity of each of the lines measured. If the continuum did not 
appear to be well-behaved, most likely because of large-scale fixed-pattern noise effects in the 
detector, then the the affected hne was not included in the fit (the SiClB MDRS measurement 
of the Ol948.686A line is such a case). Uncertainties in continuum placement are included 
explicitly in the equivalent width measurements used in the COG analysis, and implicitly in 
the minimization performed by the profile fitting. This latter approach is adequate if the 
true continuum is in fact well-represented by a low-order polynomial, and if the absorption 
being measured occupies a small fraction of the region being fit. For the broad Lyman-a 
profile a variety of continuum fits were tested to determine the sensitivity of the results to 
the adopted continuum. This will be discussed in greater detail below in the section on the 
measurement of the H I column density. 



4-2.5. Background placement errors 

The dominant background for FUSE spectra of a bright source such as HZ 43A is grating 
scattering of source fiux. The intrinsic detector background and diffuse stray light in the 
spectrograph are negligible in comparison. The grating-scattered source flux has been char- 
acterized by measuring the flux in the cores of highly-saturated H2 absorption line proflles 
in the spectra of stars with translucent clouds along the line of sight. Typical results for this 
background are that the level is about 1% of the continuum flux in the SiC channels, and 
0.3% to 0.5% of the continuum flux in the LiF channels. In the HZ43A FL^S'E' observations, 
the fiux shortward of the Lyman limit is indeed about 1% of the continuum fiux in the SiClB 
channel. None of the interstellar absorption lines in the IIZ43A spectrum are broad enough 
for the true background level to be seen in the line cores. The instrumental LSF has a broad 
component that has a FWHM of about 18 pixels, which is broad enough that even the Hi 
Lyman hne cores do not reach the level of the background. This point is discussed further 
in Section 4.2.8. The effects of uncertainty in the background level were investigated by 
performing the profile fits with the background level set to zero, set to our best estimate of 
the background, and and set to twice our best estimate of the background. The resTilting 
column density variation was typically 0.01 dex, and is included in the quoted uncertainties 
for the column densities. The background effects in the GHRS spectra are discussed in the 
section below on the Hi column density measurement. 
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4-2.6. Geocoronal emission 

Contamination of the spectrum by geocoronal emission can be significant for the large 
spectrograph apertures employed by FUSE. Such emission is quite strong at Lyman-/?, but 
declines rapidly in strength for the higher Lyman lines. Geocoronal emission is also present 
on the day side of the orbit for many of the same O I and N I transitions for which we are 
trying to measure interstellar absorption. The widths of geocoronal lines are limited by the 
size of the spectrograph apertures. For the LWRS aperture this width is typically 0.35A, or 
the equivalent of 106 km s~^ at Lyman-/3, while for the MDRS aperture this width is about 
0.045 A, or the equivalent of 13 km s~^ at Lyman-/?. 

The MDRS aperture is sufficiently narrow (~13km s~^) that the geocentric velocity 
of the interstellar gas shifts the geocoronal emission almost completely away from narrow 
lines. The red wing of the broad interstellar Lyman-/? line is filled in by geocoronal emission, 
but no difference is seen in the O I or N I absorption profiles when the day and night data 
are analyzed separately. The absorption profiles measured in the LiFl channel are actually 
deeper for the day data than the night data, presumably because the LiFl grating is more 
stable in the day than in the night and the effective resolution is correspondingly better. 

4.2.7. Uncorrected fixed-pattern noise 

Uncorrected fixed-pattern noise on small scales can lead to significant errors in the 
shapes and fiuxes of line profiles; on larger scales it can lead to continuum placement errors. 
The FUSE detectors exhibit significant fixed-pattern noise on both small and large scales, 
caused primarily by the inherent granularity of the microchannel plates (MCPs). The slowly- 
varying relative spacing of the pores in the stack of 3 MCPs gives rise to Moire fringe effects 
with a spacing of 8-9 pixels, comparable to the width of the LSF. Changes in pore spacing at 
fiber bundles and regions of dead pores can cause changes in the detector quantum efficiency 
on larger scales. The present CALFUSE pipeline flags pixels affected by dead regions of the 
detector, but otherwise fixed-pattern features are not yet corrected. 

We have mitigated the effects of fixed-pattern noise by several means. Every transition 
being examined in the FUSE bandpass is sampled by either 2 or 4 channels, so we have 
verified that each such measurement provided essentially similar results. In addition, because 
HZ 43A was observed in both LWRS and MDRS apertures, we have checked for consistency of 
the results. The MDRS data were obtained with an FP-split procedure, which substantially 
reduces the fixed-pattern noise on scales of one to two resolution elements. The resulting 
spectra typically have a S/N that is about 80-90% of that expected from photon statistics. 
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depending on the size of the region examined, primarily because the larger-scale fixed-pattern 
noise has not been removed. However, most of these fluctuations are slowly varying and can 
be removed by low-order polynomial fits to the continuum in the vicinity of absorption lines 
of interest. 



4-2.8. Uncertainties in the instrumental line spread function 

Characterization of the FUSE LSF is still in a preliminary state. There appears to be 
both a narrow and a broad component to the LSF. The narrow component is somewhat 
wavelength-dependent and is sensitive to the accuracy of the correction for grating motion 
during exposures and to the accuracy with which the user was able to co-align the spectra 
from each exposure prior to combining them. In addition, data acquired in histogram mode 
cannot be corrected for the degradation of the spectral resolution caused by rotation of 
the primary mirrors or gratings during an exposure. For example, in some channels the 
resolution is better during the night portion of the orbit and for other channels it may be 
better during the day. As a result, the effective width of the narrow component of the LSF 
will vary from one observation to another. We measured the effective width of this narrow 
component by fixing the temperature and turbulent velocity of the interstellar gas to that 
found from the GHRS measurements (see below), and letting the LSF widths vary in the 
profile fits. The results were generally consistent from one absorption line to another, and 
were typically about 9 pixels FWHM (corresponding typically to 0.062A in the LiF channels 
and O.O57A in the SiC channels). This is consistent with what is found from other FUSE 
datasets. These widths were then held fixed when performing subsequent fits. It should be 
noted that the LSF width is significantly greater near the edges of the detector active area, 
as are the uncertainties in the wavelength calibration. 

The broad component of the LSF appears to be somewhat less sensitive to the effects 
of grating rotation, primary mirror rotation, or coalignment of spectra. The initial charac- 
terization the available data indicates that about 30 percent of the area of the LSF is in 
a component with a width in the range of 17-24 pixels. Analysis of highly saturated fines 
that are flat-bottomed over regions of 150 pixels or more reveals that there is no signiflcant 
component of the LSF with a characteristic width larger than 25 pixels. The effects of this 
broad component of the LSF on the proflles of unsaturated lines may not be evident unless 
the signal-to-noise ratio of the data is very high; however, its main effect is to introduce a 
substantial apparent increase in the background under the line. If a weak line is fit with a 
single, narrow, Gaussian LSF, the background level in the fit must be increased to account 
for this effect. In the profile-fitting analysis of the HZ43A ISM absorption we have repre- 
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sented the LSF as a two-component Gaussian function. The width of the narrow component 
was determined to be about 9 pixels as described in the previous paragraph, and the broad 
component width was set to 17 pixels, with an amplitude corresponding to 30% of the total 
area. 

If the LSF has very long shallow tails that cannot be distinguished from the continuum, 
then both the profile fitting technique and curve-of-growth analysis would likely be subject 
to continuum placement errors that will cause some of the area of the line to be missed 
and therefore cause the the column density to be underestimated. However, while there is 
a component of the LSF that has roughly twice the width of the narrow component, there 
is no component with very long extended tails. In the absence of such a weak and very 
broad component, the curve-of-growth analysis is independent of the shape of the LSF, and 
thus provides a valuable check that the profile fit was not adversely affected by errors in the 
assumed LSF. 

The LSF of the GHRS modes used here is smaller than the Doppler widths of the hues 
being measured, other than for Fell, so there is no significant systematic error associated 
with uncertainties in this LSF. The Doppler width of the Fell 2382. 765A hue is comparable 
to the LSF width, so the determination of the Doppler width for this hue does have relatively 
large uncertainties. The column density, however, is not adversely affected because the line 
is so weak. 



4-2.9. Uncertainties in the f -values 

Errors in the atomic data used in the fits have the potential of biasing the results. 
If such errors are distributed randomly then the use of numerous lines will minimize any 
systematic bias of the overall column density determination. If the oscillator strengths for 
a given species are all too large or too small then the estimate for the column density will 
be biased. However, such a systematic error in the f-values seems unlikely to result in a 
consistent fit to all of the lines, especially if lines are included that arc not on the linear part 
of the curve of growth. We have tested for effects on the derived column density resulting 
from errors in the f-values by repeating the profile fits after removing each transition and 
verifying that the derived column density did not vary by more than 1 a. An examination of 
the curve-of-growth plots also shows that a satisfactory fit can be obtained that is consistent 
with each of the included points. The only transition that was found to have an f-value 
clearly inconsistent with those of the other transitions was the that of the O I 1026.473A 
line (see Sembach 2000). This line was excluded from all fits and from the curve-of-growth 
analysis. 
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4.3. Results 

4.3.1. Line of sight velocity structure 

The velocity structure of absorbing gas along the line of sight was determined from 
GHRS observations of D 1 1215.339A, Fell 2383.765A, Mgll 2796.352A, and Mgll 2803.53lA 
hne profiles. These lines are well-resolved, and are isolated from any other absorption fea- 
tures. We obtained a satisfactory fit to the lines with a single absorption component with a 
heliocentric velocity of -5.2 ± 0.1km s^^, a temperature of T = 5353 ± 948 K, and a non- 
thermal velocity of bnt = 2.0 ± 0.2km s^^ (the uncertainties are la statistical uncertainties; 
the systematic uncertainty on the absolute velocity is ± 3km s~"^). The total Doppler pa- 
rameter, b, is related to the temperature and non-thermal velocity parameter such that 
&2 — 0.01657/74 -I- The observed profiles and the best model fits are shown in Figure 1. 
There is no evidence for an additional velocity component. These values for T and bnt were 
not held fixed in subsequent fits (e.g. for O I or N l) . Instead, the derived Doppler widths for 
each species were computed and compared with the ranges allowed by the above values; in 
each case the derived Doppler widths were found to be consistent with the 1 a uncertainties 
quoted above. The b- values for O I and N I predicted from the values of T and bnt given 
above are: b(Ol) = 3.1 ± 0.2km s~^and b(Nl) = 3.2 ± 0.2km s~^ (Icr uncertainties). 



4.3.2. Hi Column Density 

The interstellar H I column density along the line of sight to HZ 43A can be determined 
by three different means: fitting the spectral shape of the EUV continuum, fitting the 
damping wings of the Hi Lyman-a profile, and fitting the converging Lyman series. This 
last method has the greatest uncertainties because the lines are typically all saturated, but 
for the low column density of this line of sight it is possible to obtain meaningful results. 

Determination of the Hi column density by analysis of the EUV spectrum of a star 
requires accurate knowledge of the stellar continuum flux. The EUV spectrum of a star is 
very sensitive to the presence of even trace quantitites of heavy elements in the photosphere, 
and substantial effort is usually required in order to model this spectrum accurately. HZ 43A 
is unusual in that no traces of heavy elements have yet been detected in its photosphere; hence 
it is not affected by this potential source of systematic error. EUVE spectra of HZ43A have 
been analyzed by several groups, using different observations, different model atmosphere 
codes, and different fitting methods. The three most recent analyses of EUVE observations 
will be described briefiy below. 
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Dupuis et al. (1995) analyzed observations from 1994 March 25,28 and 1994 May 20. 

They first determined the He I and He ll column densities by fitting the ionization edges at 
504A and 228A, respectively, and then determined the stellar T^g and interstellar H I column 
density by fits over the 400-550A bandpass. This wavelength range was chosen because it 
was least sensitive to variations in surface gravity or heavy clement abundance. The fitting 
was done with model spectra computed from pure H model atmospheres (Vennes & Fontaine 
1992). The results were T^s = 51100 K ± 500 K, and Nhi = (8.7 ± 0.6) x 10^^ cm'^ (la). 
They also estimated that systematic uncertainties in the effective area of 20% would change 
Teff by 3% and Nhi by 10%. Uncertainties in the V magnitude (used to normahze the 
models) would have similar effects, but are small in comparison. The ionization fraction of 
hydrogen was estimated by assuming that the ratio of total H to total He was 10:1, that 
there was no appreciable He III, and comparing the measured Nh i to the measured He I and 
Hell column densities. The resulting ionization fraction for H, Jh = N(Hll) / N(H), was < 
0.4 ±0.1. 

Barstow et al. (1997) analyzed the 1994 May 20 observation also, but used H+He 
models computed by Kocstcr (1991) and determined all stellar and ISM parameters from a 
simultaneous fit to the entire EUVE spectrum. The values for Toff and log g were constrained 
to remain within the 1 a uncertainties of the values determined by Napiwotzki et al. (1993). 
The V magnitude used, 12.99, corresponds to model fluxes that are fainter by a factor of 
1.077 than would be the case if the present best value of V = 12.909 were used (Bohhn 
2000) . This would result in an underestimate of the derived Nh i column density by a few 
percent. They performed flts with two classes of models: homogeneous H+He mixtures, and 
stratifled compositions. The derived interstellar H I column densities and 1 a uncertainties 
were Nhi = (8.8 ±0.2) x lO^^cm'^ (homogeneous) and Nhi = ^-^t^l x lO^'^cm^^ (stratified). 
The difference between the results for the two models exceeds the statistical errors in the 
fits, indicating that systematic uncertainties in the nature of the appropriate model are 
significant. The ionization fraction was derived in the same manner as by Dupuis et al. 
(1995) to be fn = O-lQto.ll for homogeneous models, and fn = 0.22l[J:i4 for stratifled 
models. We have assumed that the results presented by Barstow et al. (1997) supercede the 
earher results in Barstow et al. (1995), so we will not discuss the earlier results here. 

Wolff et al. (1999) analyzed observations obtained 1997 June 25,26. They flrst flt the He I 
and He ll ionization edges to determine the corresponding column densities, and then flt the 
full spectrum to determine the H I column density. The stellar spectrum was computed using 
the LTE model atmosphere code of Koester (1996), assuming a pure hydrogen composition. 

The atmosphere parameters were not free, but fixed at values derived from flts to the Balmer 
lines. The values used for Tcfi were 50800 ± 300 K, and the corresponding column densities 
were Nhi — 8.8511q;3qo xlO^^ cm~^. The uncertainties in the column densities are dominated 
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by the uncertainty in the effective temperature (Wolff 2001, private communication). The V 

magnitude used to normaUze the models was 12.89; this would overestimate the stellar flux 
by about 2% and would result in only a very slight overestimate of the Hi column density. 
The ionization fraction was derived in the same manner as the previous two analyses to be 

fn = 0.03lo:o3- 

An EUV spectrum of this star was also obtained by HUT and analyzed by Kimble 
et al. (1993) to derive Nhi = (6.5 ±0.5) x lO^W'^ (la statistical errors; systematic 
uncertainties were estimated at 15-20% and are dominated by calibration uncertainties). 
This falls below the more recent EUVE determinations by approximately the combined 
random and systematic uncertainties. These uncertainties are considerably larger than those 
of the much higher signal-to-noise EUVE spectra, so we will not consider this measurement 
further. The three EUVE measurements described above will be addressed again at the end 
of this section. 

We fit the GHRS spectrum of the interstellar Hi plus Dl Lyman-a profile towards 
HZ 43A to determine the properties of the neutral hydrogen in this direction. The spectrum 
was first scaled by a factor of 0.866 in order to match our model fiux over the regions 
1212.25 - 1213.25A and 1218.0 - 1218.5A. The spectrum was then normalized by an NLTE 
stellar model prior to fitting (see Section 4.1 for a description of the model). We found that 
additional normalization by a low-order polynomial was required to account for potential 
small wavelength-dependent errors in the flux calibration. The apparent photospheric radial 
velocity of HZ43A {vph) was measured by Reid (1996) to be +21 km s^^. The spectrum and 
the model (shifted by +21 km s^^) are shown plotted in Figure 2. A visual inspection of the 
figure shows that at least a linear polynomial will be required as an additional correction to 
the fiux calibration. 

No explicit uncertainty was given by Reid (1996) for the radial velocity, but it is probably 
about 10-15 km s^^. We therefore examined the effects on the derived Hi column density of 
varying this velocity. We created stellar models with the photospheric velocity set at zero, 
at the nominal value of +21km s""^, and at ±5, ±10, ±15 km s~^ with respect to the nominal 
value. The GHRS Lyman-a spectrum was normalized by each model in turn, and the best 
fit to the spectrum was generated for fixed values of log Nm ranging from 17.80 to 18.00 
in steps of 0.01 dex. The resulting contours are plotted in Figure 3. The Fell and Mgll 
lines were fit simultaneously with the Lyman-a H I and D I line profiles. Because the Lyman- 
a proflle extends over several Angstroms, it may be affected more by uncertainties in the 
flux calibration than are narrow lines. To investigate this possibility, we generated these 
contours for continuum fits performed with both fourth-order and seventh-order polynomials. 
There was little correlation in either case between log Nhi and the photospheric velocity. 
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For the fourth-order fits, the minimum was obtained for a photospheric velocity of 26 km 
s~^, but the increased slowly for smaller velocities and increased rather steeply for higher 
velocities. For the seventh-order continuum fits, the minimum was also at 26 km s~^, but 
it increased rather rapidly for both lower and higher velocities. The H I column density and 
2 a uncertainty derived was log Nhi = 17.895]'^o'o35 from the fourth-order continuum fits, 
and log Nhi = 17.905l;Q;Qg5 from the seventh-order fits. The corresponding values were 
2750 for 2871 degrees of freedom, and 2711 for 2868 degrees of freedom, respectively. Based 
on these fits, we have adopted the value log Nhi — 17-901° Qg for our combined GHRS and 
FUSE datasets. 

We examined the behavior of the continuum fits by varying the order of the polynomial 
from 1 to 8, for photospheric velocities {vph) of 0, -|-5, -1-11, +21, and -|-31km s~^, and 
allowing the Hi column density to vary freely in each fit. For velocities of and -|-5 km 
s"^ there was relatively little effect of polynomial order on the reduced (for orders > 
1). For higher velocities the x^ improved significantly for orders up to 4 and improved only 
slightly thereafter (with minima at order 7), except that for -|-31km s~^ the x^ dropped 
again for order 8. An examination of the polynomials in each case show that the dominant 
term gradually changes from linear at low Vph to cubic at high Vph, and that the amplitudes 
increase from ± 2% at low Vph to ± 3% at higher Vph. The higher-order terms gave slight 
improvements to the overall x^, but did not change the character of the continuum fits. 

We evaluated the effects of the uncertainties of the background flux level on the derived 
log Nh I- Our best estimate of the background level at Lyman-a was 0.5% of the continuum 
flux level, obtained by letting the background be a free parameter in the flt. We then varied 
the background level by +2% and -2% of the continuum flux level, causing log Nh i to change 
by only ± 0.012. 

We were unable to obtain a satisfactory flt to the observed Lyman-a profllc with a single 
component of neutral material, but found that a good flt was possible only if a second, low- 
column density absorber were included. The temperature, column density, and relative 
velocity of this component were free parameters in the fits to the Lyman-a profile used 
to construct the x^ contours described above. For the fits in which the continuum was 
represented by fourth-order polynomials, the column density of this second component was 
log Nhi = 14.91q5 (2 cr)- For the fits in which the continuum was represented by seventh- 
order polynomials, the column density of this second component was log Nhi = 14.9lg2 (20"); 
We will adopt the value log Nh i = 14.91q'5 (2 '^)- The temperature of this component was 
30,000 K ± 10,000 K (2 a) for both cases. As the column density of the main Hi component 
was increased, the derived column density for this second component would also increase, 
and its temperature would decrease. Similarly, as Vph was increased, shifting the NLTE 
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core of the photospheric Lyman line closer to the red side of the saturated portion of the 
interstellar absorption profile, the column density of this second component would increase 
and the temperature would decrease. The heliocentric velocity of this component was found 
to be about -2.0km s~^, but as Vph was varied from it would shift in the opposite direction 
(i.e. it would increase as Vph decreased). We suggest that this low-column density, high- 
temperature Hi component may be due to the "hydrogen wall" about the Solar System. 
The hydrogen wall is an accumulation of gas at the boundary of the heliosphere, where the 
solar wind and LISM interact. Hydrogen walls have been detected around the Sun along 
several lines of sight, and around several other stars (Wood & Linsky 1998; Linsky & Wood 
1996; Wood et al. 2000). The gas in these walls is fairly hot, 30,000K - 180,000K, and Hi 
column densities range from log Nhi — 13.75 to log Nhi = 14.75. The properties of the 
second Hi component in our fits are consistent with these other detections. One potential 
problem with this interpretation, however, is that the Lyman-a profile observed for the 
similar high-latitude sightline to 31 Com (/ = 115°, 6 = +89°) can be fit very well without 
any high-temperature component (Wood et al. 2000). The models of the interaction of the 
hehosphere and the LIC described by Wood et al. (2000) do predict absorption by such a 
hot component for the line of sight to 31 Com, so one possible explanation is that the LIC 
is inhomogeneous on small scales and that the corresponding conditions at the boundary of 
the heliosphere vary accordingly. 

The best fit to the normalized Lyman-a profile, including the main component and 
the high-temperature component, is shown in Figure 4. The relative velocity of the main 
absorbing component matched that found for the Dl, Mgll, and Fell lines, about -5km 
This differs significantly from that expected for the local interstellar cloud (LIC), -8.8km s~^, 
and the column density is much higher than would be expected from the LIC along this line 
of sight (log Nhi = 16.47; Redficld & Linsky (2000)). The velocity differs even more from 
that expected for the 'C cloud, -12.1km s~^. The velocity is consistent, however, with that 
found for the single absorber along the line of sight to 31 Com by Piskunov et al. (1997) (-4 
± 1km s~^). This absorber is called the North Galactic Pole cloud by Linsky et al. (2000, 
see their Figure 1). 

The FUSE spectrum of the converging Lyman series might in principle be used to try to 
determine the H I column. Such fits give results consistent with the column density derived 
from the Lyman-a profile, but with rather large uncertainties. The main difficulty is that the 
fiux calibration is inherently uncertain at precisely the wavelengths needed to perform the 
fits to the Lyman edge, because there is no reference star available for which unattenuated 
continuum fiux is present at the Lyman edge. For any given set of parameters we find that 
the fiux is overestimated at some wavelengths and underestimated at others, by about ±5%. 
If we allow for such uncertainties we find that variations of the gas Doppler b parameter 
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within 1 a limits result in variations in the derived log Nhi of about 0.15 dex, which is 
too large to provide a significant constraint on our results. Ultimately, we may be able to 
determine a more accurate calibration at these wavelengths by means of combining data from 
several white dwarfs for which reliable Hi columns have been determined by other means. 
The interstellar Hi absorption at the Lyman edge predicted from the column densities, 
temperatures, and turbulent velocity derived above from the GHRS data is shown in Figure 
5. The agreement is excellent, indicating in particular that the Dopplcr width of the gas is 
accurately known. The effects of the broad component of the instrument LSF is clear in this 
Figure, as the flux in the cores of the saturated Lyman lines is several times higher than the 
background level seen below the Lyman edge. 

We will combine the available £'C/VE-derived measurements of the Hi column density 
with our GHRS / FUSE measurement to provide our best estimate of Nhi. The EUVE 
measurement uncertainties are driven in large part by the uncertainty in the effective tem- 
perature. Each group used different estimates for this uncertainty, but they show similar 
sensitivity in the sense that a change in Tcfj of 500K would correspond roughly to a change 
in Nhi of 0.6 x 10^''cm~^ (this scaling is not apparent in the Barstow ct al. 1997 paper, 
but in the earlier Barstow et al. 1995 paper Nhi appears to be roughly 50% more sensitive 
to changes in Tgff). If we take the uncertainty for Teg from Dupuis et al. (1995) as rep- 
resentative, the corresponding uncertainty in Nhi from an individual EUVE measurement 
is roughly equal to the 1 a uncertainty in our determination of Nhi from the GHRS and 
FUSE data. Therefore, we will simply take the mean of the four measurements (after flrst 
averaging the two values obtained by Barstow et al. (1997) for both the homogeneous and 
stratified models), which gives Nhi = 8.511 x lO^^cm"^, or log Nhi = 17.930. The standard 
deviation of the four values was 0.398 x 10^^, or 0.02 dex, but this may not be an adequate 
representation of the true uncertainty. Each of the EUV measurements is affected in a simi- 
lar fashion by the uncertainties in the effective temperature. Barstow ct al. (1997) were able 
to obtain relatively small uncertainties for Nh i in the context of either the homogeneous or 
stratifled models, but the differences between the two corresponding column densities was 
0.5 X lO^^cm"^, or almost as much as for the nominal uncertainty of 500 K for Tes. It is 
worth noting that the two recent optical determinations of Teff have uncertainties greater 
than 500 K: 49000 ± 2000 K {3 a) by Napiwotzki et al. (1993) and 50822 ± 639 K (1 o") by 
Finley et al. (1997). Similarly, as we showed above in Section 2 the standard deviation for 
all the recent determinations of Tcs is about 700 K. Therefore, we prefer to use a somewhat 
more conservative estimate for the uncertainty in Nhi, and adopt the value log Nhi ~ 
17.930 ± 0.060 {2a). 
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4-3.3. Ol Column Density 

The FUSE bandpass contains numerous O I absorption lines of varying strength. The 
total column density to HZ 43A is sufficiently low that even the strongest O I lines are close 
to being on the linear part of the curve of growth. The data were analyzed with profile 
fits and a curve-of-growth analysis. Good quality data were obtained in all channels for 
the MDRS observation. Separate analyses of the day and night MDRS data actually give 
a somewhat higher value for N(Ol) than the night data; the opposite would be expected if 
airglow contamination were significant. 

The MDRS data were analyzed in 4 ways: the program Owens.f was used by two mem- 
bers of our group to perform independent profile fits, and a third member of our group 
performed both a profile fit using the code of Fitzpatrick (Spitzer & Fitzpatrick 1993) and 

a curve-of-growth analysis. Each analysis was performed with different choices of contin- 
uum regions to fit, different choices of polynomial order for the continuum fits, and different 
combinations of background level and LSF. The resulting column densities, with 2 a uncer- 
tainties, were: log Nqi = U.A7to^l, log ^oi = U.AStofj, log ^oi = 14.50 ± 0.04, and 
log Nqi = 14.5llQ;Qg. The corresponding b- values and 2 a uncertainties from each fit were: 
2.8 ± 0.2, 3, 2.5 ± 0.4, and 2.7±0.3. The b-value for the second fit is poorly-determined 
because only unsaturated lines were included, which are insensitive to b. The first two b- 
values were from simultaneous fits to Ol, Nl, and Nil lines in the FUSE data; if only the 
O I lines are included in the first fit the uncertainty increases to ±0.40 km s""^. The b- values 
fall somewhat below the value predicted from the line of sight velocity structure derived 
from the GHRS data, but are consistent within 2 a. The above determinations are in good 
agreement with one another; we adopt the value log Nqi = 14.49 ± 0.08 {2 a). A sample of 
the fits to the O I lines is shown in Figure 6, and the curve of growth is shown in Figure 7. 

A number of tests were performed to examine the magnitude of the various possible 
systematic effects described in Section 4.2. Increasing the background level by a factor of 
2 in the SiC channels and a factor of 4 in the LiF channels increased the value of log Nqi 
by 0.02 dex, decreasing the background by a factor of 2 decreased log Nqi by < 0.01 dex. 
Fitting only detector 1 data or only detector 2 data, to test for fixed-pattern noise effects 
and detector-dependent resolution or background effects, changed log Nqi by < 0.01 dex. 
Eliminating any single line from the fit, to test for erroneous f- values or fixed-pattern noise 
effects, had no effect on log Nqi- Ehminating both the fine at 1039A and the multiplet at 
988A from the first profile fit causes log Nqi to increase by 0.04 dex, but the uncertainty 
of the result doubles. Eliminating these lines from the COG analysis causes the resulting 
log Nqi to increase by only 0.02 dex. The LSF was varied by doubling the amplitude 
of the broad component, which increased log Nqi by 0.02 dex, and eliminating the broad 
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component (with a corresponding increase in the background level to match the cores of the 
saturated Lyman lines), which decreased log Nqi by 0.02 dex. These effects have all been 
included in the uncertainties quoted above for the first profile fit. 

4.. 3. 4.- Nl Column Density 

As with Ol, the FUSE bandpass contains numerous Nl absorption lines of varying 
strength. In this case the low column density to HZ43A means that many of the hues 
that are ordinarily the most useful arc not detected. High quality data were obtained in 
all channels in the MDRS observation, and usable data were obtained in all channels in the 
LWRS observation. We analyzed both the complete dataset and the pure night data for the 
LWRS observation and found that the gcocoronal emission from N I during the day portion 
of the orbit did not significantly affect the absorption line profiles. The lines strong enough 
to be useful in the FUSE data were the triplet at 1134A and the lines at 953.655A and 
953.970A. 

The MDRS data were analyzed by having two members of the group perform indepen- 
dent profile fits and a third member perform a curve-of-growth fit. The curvc-of-growth anal- 
ysis also included the Nl multiplet 1199.550A, 1200.223A, and 1200.710A from the GHRS 
G160M dataset. For the MDRS data, the results were log Nni = 13.51 1^;^^, log A^jvj = 
13.5llo!o5) log Nni — 13.511q;q5 (2 a uncertainties). The b-value and 2 a uncertainty 
derived from the curve of growth analysis was S.Olo rkni s~^; no determination of b was 
possible from the profile fits because the lines used were all on the linear part of the curve 
of growth. The adopted value for log i is 13.51 ± 0.06. The fit to the MDRS Nl profiles 
is shown in Figure 8, and the curve of growth is shown in Figure 9. 

The Nl fits were subjected to the same tests for systematic errors as described in the 
previous section for O I. Because all of the N I lines in the FUSE band relatively weak, they 
are less sensitive to variations in the LSF or the background than partially saturated lines; 
none of these tests caused a variation in the log iVjv/ greater than 0.01 dex. 

4.3.5. Dl Column Density 

Because of the low Hi column density along this line of sight, absorption by D I is 
expected to be detectable at Lyman-a, Lyman-/?, and possibly at Lyman-7. The D I Lyman- 
a profile is shown in Figures 2 and 4; it is clearly well-resolved. The D I Lyman-/? absorption 
is clearly seen in all FUSE channels in the MDRS observation, but is clearly detected only in 
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the night portion of the LWRS observation in the LiFl channel. The MDRS SiC2 Lyman-/3 
spectrum was rather noisy, however, and was not included in the fits. D I Lyman-7 absorption 
is not detected in either SiC channel. 

The GHRS Lyman-o; profile was first fit independent of the FUSE data. The resulting 
value for the column density was log Ndi — 13.16 ± 0.04 {2a), with a b- value of 7.1 (implying 
the line is wcll-rcsolvcd by the GHRS echelle resolution). Changing the background level 
from our best estimate by ± 2% of the continuum level changed log N£,i by ±0.02 dcx. 
A simultaneous fit to the GHRS Lyman-a spectrum and the FUSE MDRS Lyman-/? data 
from channels LiFl, LiF2, and SiCl gave log Ndi = 13.15 ^0945 (2o"). Varying the fitting 
parameters for the FUSE data as described in Section 4.3.3 above for O i resulted in variations 
in log Ndi of ±0.01. Profile fits to the FUSE MDRS data alone gave log Ndi = 13.11 ± 
0.20. The LWRS LiFl Lyman-/3 spectrum was not included in the fits. The Dl line is 
detected in this spectrum, and is consistent with the fits to the other data, but the uncertain 
corrections for airglow contamination in the LWRS spectrum, even for the night only data, 
means that these data will not significantly constrain the overall fit. The FUSE D I line 
profiles are shown in Figure 10, along with the best fit to the combined FUSE and GHRS 
datasets. 



4-3.6. Other Species 

A number of other species were detected in the ISM along the line of sight to HZ 43A. 
These include Cll, Cm, Nil, Ovi, Sill, An, and possibly Nlll in addition to the Mgll 
and Fell mentioned earlier. The column densities for these species are listed in Table 5. 
Of particular interest is the column density for Nil (log A^jv//= 13.62]|;q}7), which slightly 
exceeds that of N I. The ionization potential for N I is greater than that for H I, so ordinarily 
one expects N to be predominantly neutral. The local ISM is dominated by hot gas that may 
be far from ionization equilibrium, however, so estimates of the total abundance of nitrogen 
in the local ISM must account for possible large ionization corrections to the measured 
Nl abundances. The Nil transitions at 915.613A and 1083. 993A are ordinarily strongly 
saturated, but the total column density to HZ43A is low enough that an accurate Nil 
column density can be determined. The Nil absorption lines and best-fit line profiles are 
shown in Figure 11. The NlII column density quoted in Table 5 is an upper limit because 
the NlII 989.799A is blended with the Sill 989.873A line. If the Sill column density derived 
from the 1193. 290A line and the b-value derived from Ol are used to predict the equivalent 
width of the Sill 989.873A line, then it is likely that at least half, and perhaps all, of 
the measured equivalent width at 989. SA is caused by Sill rather than Nlll. Even if the 



-25- 



measured absorption results primarily from Nlll, the upper limit suggests that Nlll is a 
minor contributor to the total N column along this sightline. 

No H2 was detected. The range of possible column densities was determined by using 
the profile-fitting program Owens to examine the regions of the spectrum in the vicinity of 
the two-dozen strongest H2 absorption lines that were not blended with any other features. 
The column density of H2 was then gradually increased until the change in chi-squarcd was 
unacceptable. The resulting 3 — a upper limit to the column density of H2 was 1.2 x 10^^ 

_2 

cm . 

The column density for C ll shown in Table 5 is larger than that for O I, possibly suggest- 
ing that C is over-abundant on this sightline, that O is depleted, or that there is substantial 
component of Hll gas along this line of sight. The large column density of N II supports the 
latter explanation. It should be noted also that the column density for ClI is particularly 
uncertain, because the saturation effects are large and depend on the assumption that the 
b- value for C II is the same as that of O I. If there is a substantial contribution to the column 
density from H ll gas, the effective b- value for C ll will likely be larger than that for O I and 
the C II column density will be overestimated accordingly. This source of systematic error 
is not included in the uncertainties quoted in Table 5; the column density derived from the 
apparent optical depth method (log N(Cll) = 13.63) is a reasonable lower limit. 

We note that Cll* at 1037.018 A is not detected in our FUSE observations. The 3a 
upper limit to the equivalent width of this hne in the LiFlA data is Wx < 5.5 mA, including 
a crude estimate of the continuum placement uncertainties. This imphes log A^(Cll*) < 12.7 
{3a) along the HZ 43A sight line. If one assumes the excitation of C ll* is due to electrons, 
the electron density may be calculated by comparing the column densities of the excited- 
and ground-state C ll: 

ivrcii*) 

(see Spitzcr & Fitzpatrick 1993). Assuming T = 5353 K (sec Section 4.3.1) and that 
logiV(Cll) = 14.83 (derived using the Ol 6-value), we find ng < 0.1 cm~^. Given the 
uncertainties in the latter assumption, however, a more conservative value is ng <^ 1.5 cm~^, 
derived using the lower limit for the C ll column density determined with the apparent optical 
depth method. 
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5. Summciry 

The line of sight to HZ43A was investigated using a combination of FUSE spectra, 
archival GHRS spectra, and results obtained from EUVE spectra taken from the literature. 
The line of sight to HZ 43A is particularly simple, with only a single velocity component 
discernible. An additional high-temperature component was found with a very low column 
density, detectable only in H I. The hne of sight to HZ 43A exits the LIC after a very short 
distance, and the velocity of the main component of the absorbing gas is not compatible with 
the velocity of the nearby G cloud in the LISM, hence this material is presumed to be located 
in the North Galactic Pole cloud. Doppler widths for each species determined from profile 
fits or curve of growth analyses of the FUSE data were consistent with the temperature and 
turbulent velocity for the gas derived from the GHRS high resolution data. Consistent fits to 
the D I column density were obtained when fitting the FUSE and GHRS data separately and 
simultaneously, indicating that the Dl Lyman-a hne is not affected by saturation effects. 
The FUSE spectra were used to determine O I and D I column densities by both profile fitting 
and curve of growth analyses, with consistent results obtained in each case. Careful attention 
was paid to numerous potential sources of systematic error. Our adopted H I column density, 
log Nh 1= 17.93 ± 0.06 (2 cr), is the mean of our measurement from the GHRS Lyman-a 
profile and the EUV measurements of Dupuis et al. (1995), Barstow et al. (1997), and Wolff 
et al. (1999). 

Our results are summarized in Table 6. All uncertainties in the table are 2 a (but 
values quoted below from the literature are 1 a). Our result for D/H along the line of 
sight to HZ43A is (1.66 ± 0.28) x 10~^. The GHRS dataset analyzed in this work has been 
previously analyzed by Landsman et al. (1996), who obtained the similar value for D/H of 
1.6 X 10^^. Our value for D/H along this sightline exceeds by approximately 1.5 ex both the 
mean value reported by Linsky (1998) for measurements of D/H within 100 pc of the sun, 
(1.47 ± 0.10) X 10~^, and the mean value of the FUSE measurements reported in this first 
set of papers (Moos et al. 2001). 

Our result for Ol/Hl is (3.63 ± 0.84) x 10~^. Converting this ratio to a logarithmic 
abundance log (O l/H l) + 12.00 gives 8.56lo;ii- The ionization balance for both O and H are 
linked by resonant charge exchange reactions (see Jenkins et al. 2000 and references therein), 
and there is no apparent H2 present, so Ol/Hl should be representative of 0/H along this 
line of sight. Meyer et al. (1998) report an average gas-phase abundance for 13 sight fines 
of (3.43 ± 0.15) X 10~^ (after correction for our preferred Ol A1355 oscillator strength of 
/ = 1.16 X 10^^ from Welty et al. 1999), and they estimate no more than 0/H~ 1.8 x 10""^ 
can be incorporated into grains, implying an upper limit to the total (gas+dust) ISM oxygen 
abundance of 5.2 x lO^''. Our 0/H measurement is consistent with the average gas-phase 
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value of Meyer et al. (1998). 

Sofia & Meyer (200 lb, a) have demonstrated that the value of 0/H determined from F 
and G star photospheres, (4.45 ± 1.56) x 10~^, is consistent with two recent determinations 
of the solar abundance, (5.45 ±1.00) x 10"^ (Holweger 2001) and (4.90tg|^) x 10"^ (AUende 
Prieto et al. 2001), and with the total oxygen abundance estimated from the work of Meyer 
et al. (1998). Sofia & Meyer (2001b, a) argue that the recent solar system values should there- 
fore be used as the ISM abundance standard for oxygen. Our value for 0/H in the gas phase 
along the line of sight to HZ43A is significantly lower than each of these values, although 
not by more than 2 a (given the large uncertainties in the solar system determinations). Our 
result is in agreement with the compilation of B-star abundances presented by Sofia & Meyer 
(2001b, a), (3.50 ± 1.33) x 10~^, although these authors argue that B-star abundances are 
unhkely to be appropriate measures of the ISM standard. It seems hkely that the gas-phase 
oxygen abundance measured by Meyer et al. (1998) implies modest depletion of oxygen in 
the ISM towards more distant stars; given the agreement of our 0/H measurement with the 
Meyer et al. average, it seems likely a similar degree of oxygen depletion is present in the 
cloud observed along the HZ 43A sight line. 

The ionization fractions for N and H are not as strongly coupled as those of O and H 
(Jenkins et al. 2000), so it is possible for N to be more highly ionized than H, and therefore 
for the gas-phase abundance of N relative to its solar system value to be less than that of 
O. If we combine our column densities for Nl, Nil, and the maximum likely column density 
for NlII (assuming only half of the measured absorption is due to Sill), we obtain a total 
nitrogen column density of = 7.88^^:43 x 10^^ and a ratio N/Hl = (9.26 ± 1.87) x 10-'\ 
This ratio provides an upper limit to N/H, corresponding to the limit that that hydrogen 
along this line of sight is predominantly neutral. Given the significant column density of 
N II, it is reasonable to expect that there will also be a significant column density of H II 
present. Converting to a logarithmic scale with log N^i = 12.00 gives an abundance for N 
of 7.97lo5'o- This is consistent with the solar abundance value of Grevesse & Noels (1993), 
7.97, and the more recent value from Holweger (2001) of 7.93 ± 0.11, and somewhat exceeds 
the value of 7.81 for nearby B stars (Gies & Lambert 1992). We can set a lower limit to 
the N abundance by setting the ionization fraction of H to the upper limit of 50% found 
by Dupuis et al. (1995) (the other EUVE analyses obtained lower ionization fractions for 
hydrogen). The resulting logarithmic abundance is 7.79lo;iO) essentially equal to the value 
for B stars, and consistent with the ISM gas phase abundance of Meyer et al. (1997). Thus 
N is not significantly depleted onto dust grains along this line of sight. 

Further discussion of the implications of these results and comparisons with other sight- 
fines can be found in the companion paper by Moos et al. (2001). 
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Fig. 1. — Absorption profiles of Dl, Fell, and Mgll are plotted as a function of heliocentric 
velocity. The histograms are the GHRS echelle measurements, and the smooth line is the 
single-component fit. No evidence for a second velocity component is seen. The D I line sits 
on the broad wing of the interstellar H I Lyman-a profile; the truncated peak at +65 km s~^ 
in the D I plot is geo coronal Lyman-a emission. 
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Fig. 2. — The GHRS echelle spectrum of the Lyman-a profile is shown as the histogram, 
the stellar model as the smooth curve. The model is shifted to a heliocentric velocity of 
+21 km s^^. The positions of the interstellar Hi and Dl absorption lines are marked. The 
truncated peak at 1215. 6A marked with an earth symbol is geocoronal H I Lyman-a emission. 
The GHRS spectrum was normalized to match the mean flux of the model over the regions 
I212.25A-I213.25A and I2I8.0A-I2I8.5A; additional low-order polynomial corrections were 
applied to the continuum flux when fltting (see text). The zero level is shown as the dashed 
line. 
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Fig. 3. — Plots of the contours are shown for fits to the GHRS Lyman-a profile as a 
function of log Nh i and photospheric velocity. The top plot shows the results for 4*^^ order 
polynomial continuum fits, and the bottom plot shows the results for 7*^^ order polynomial 
continuum fits. 
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Fig. 4. — The normalized GHRS spectrum of the Lyman-a profile is shown as the histogram, 
and the best fit is shown as the smooth curve. The thin curves show the continuum fit and 
the individual Hi and Dl profiles; the dashed curve shows the low-column density high- 
temperature Hi component; and the thick solid line shows the overall fit. The stellar model 
in this fit was shifted to a heliocentric velocity of 26 km s~^, which provided slightly better fits 
than the nominal value of 21 km s^^. The positions of the interstellar Hi and D I absorption 
lines are marked. The truncated peak at 1215. 6A marked with an earth symbol is geocoronal 
H I Lyman-a emission. The zero level is shown as a dashed line. 
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Fig. 5. — The FUSE MDRS SiClb spectrum of the converging Hi Lyman series is plotted 
as the histogrammed points, the profiles resulting from the best fit to the Lyman-a profile 
as the smooth line. The assumed continuum is shown as the thin solid line. The LWRS 
spectrum is essentially identical, but with lower signal to noise. 
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Fig. 6. — The observed FUSE MDRS Ol line profiles are shown as the histogrammed data 
points, the continuum fits are the thin solid lines, and the profiles computed from the adopted 
O I column density and doppler width is plotted as the thick solid curve. The data are un- 
binned. 
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Fig. 7. — The equivalent widths of O I hnes measured from FUSE MDRS data are plotted, 
along with the best-fit curve of growth (b = 2.72 km s~^). The inset shows contours of 
constant Ax^ in the b-log N plane. The measurements of lines from multiple channels have 
been offset slightly in log NAf to avoid confusion. The wavelengths of the lines are marked 
along the top of the figure. The plotted errorbars indicate the la measurement uncertainties 
in the equivalent widths. 
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Fig. 8. — The FUSE MDRS N l line profiles are shown as the thin histogram and the best 
fit is shown as the thick sohd line. The continuum fit is the thin smooth fine. The SiC data 
have been binned by two pixels for plotting purposes. 
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Fig. 9. — As for Figure 7, but for Nl. The solid circles are from FUSE MDRS data, the 
open circles from intermediate-resolution GHRS spectra. 
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Fig. 10. — The FUSE spectra of Hi and Dl Lyman-/3 and Lyman-7 absorption are plotted 
as the thin histogram, and the best fit to the combined GHRS and FUSE data as the thick 
smooth fine. No D I absorption is seen at Lyman-7, but it is clearly seen at Lyman-/3. The 
geocoronal Lyman-/9 emission is clearly seen in the red wing of the interstellar absorption; it 
is well-separated from Dl. Geocoronal emission at Lyman-7 is present, but quite weak. The 
Ol 1026. 473A absorption predicted from the measured Ol column density is also plotted; 
no absorption is seen at this wavelength so it has been assumed that the oscillator strength 
for this line is in error. It was not included in any fits. 
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Fig. 11. — The FUSE MDRS N ll line profiles are shown as the thin histogram and the best fit 
is shown as the thick solid line. The continuum fit is the thin smooth line. The neighboring 
interstellar Hi line profiles are labelled with tick marks in the plot of the 915.663Aline. 
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Table 1. Target Summary for HZ 43A 



Quantity 


Value 


Reference 


Spectral Type 


DAI 


1 


0^2000'^ 


13''16™21!78 


2 


A ^ 
"2000 


+29°05'55'.'5 


2 


I 


54.11° 


2 


b 


+84.16° 


2 


(pc) 


68 ±13 


3 


V 


12.914 


2 


B-V 


-0.31 


2 


TeS (K) 


50,900 


1 


log g (cm s~^) 


8.0 


1 



^Corrected for proper motion to epoch 
2000.5 

'^Trigonometric parallax 

References. — (1) Finlcy et al. 1997; (2) 
Bohlin et al. 1995; (3)van Altena et al. 1995. 
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Table 2. Log of FUSE observations. 



Dataset ID 


Aperture 


Obs date 


Texp (total)'' 


Texp (night)" 


N ^ 
exp 


M1010501 


LWRS 


Feb 19 2000 


6092 


2333 


8 


P1042301 


LWRS 


Apr 22 2000 


14447 


5136 


23 


P 1042302 


MORS 


Feb 08 2001 


39606 


20572 


85 



''Exposure duration in seconds. 
''Number of individual exposures. 
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Table 3. Log of GHRS observations. 



Dataset 


Grating 


Aperture and 


Wavelength range 


Start 


Exposure 


ID 


and order 


substep pattern 


A 


Time (UT)^ 


Time (s) 


Z2R50104T 


Ech A - 


46 


LSA 7 


1212-1218 


00:21 


7834 


Z2R50107T 


Ech A - 


46 


LSA 7 


1212-1218 


05:05 


7181 


Z2R50204T 


Ech B - 


24 


LSA 7 


2376-2387 


09:59 


2176 


Z2R50205T 


Ech B - 


20 


LSA 7 


2791-2807 


11:27 


2176 


Z2R50209T 


G160M 


SSA 5 


1192-1228 


13:10 


1414 



''AH exposures began on July 30 1995. 
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Table 4. Equivalent Width Measurements 



Species A/ logA/^ Channel Wx [mA]^ S/N^^ 



FUSE 



C II 1036.337 2.088 

C III 977.020 2.870 

N I 953.655 1.376 

N I 953.970 1.521 

N I 1134.165 1.238 

N I 1134.415 1.528 

N I 1134.980 1.693 

N II 915.613 2.180 

N II 1083.994 2.097 

N III 989.799 2.085 

O I 936.630 0.534 

O I 948.685 0.778 

O I 971.738 1.123^ 



LiFlA 


37. 


0+1.7 


42 


LiF2B 


39. 


8+2-2 
■0-2.8 


26 


SiClA 


39. 


4+2-2 
■^-2.4 


19 



SiClB 14.6l|^ 16 
SiC2A 15.61};^ 22 



SiClB 


4.6lf^ 


15 


SiC2A 


4.711? 


17 


SiClB 


9.611^ 


16 


SiC2A 


7 9+1.9 


20 


LiFlB 




47 


LiF2A 


4+1.0 

^•+-0.9 


47 


LiFlB 


8 9+^-'^ 


46 


LiF2A 




43 


LiFlB 




46 


LiF2A 


11.11- 


43 



SiClB 20.31^:1 13 
SiC2B 23.311? 14 

SiC2A < 8.2l2;?e 25 



SiClB 


8 Q+2-9 
o.O 2 7 


13 


SiC2A 


11.4+1? 


15 


SiClB 


7.81^7 


12 


SiC2A 


8.8li| 


15 


SiClB 


15.91^:^ 


10 
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Table 4 — Continued 



Species A/ logA/^ Channel Wx [mA]" S/N'^ 

O I 976.448 0.509 

O I 988.655 0.914 

O I 988.773 1.662 

O I 1039.230 0.980 



SiC2A 


16.511^ 


13 


SiClB 


8.4t|I 


15 


SiC2A 


8.0i;° 


18 


SiClB 




20 


SiC2A 


is.otll 


23 


LiF2B 


15.4t|^ 


14 


SiClB 


27.5t|^ 


20 


SiC2A 


30.41^:^ 


24 


LiF2B 


26.ltli 


14 


LiFlA 


19.1+^2' 


42 


LiF2B 


16.9t|^ 


26 


SiClA 


16.411^ 


19 



VI 


1031.926 


2.137 


LiFlA 
LiF2B 


8.51- 
8.41- 


39 
25 


Ar I 


1048.220 


2.442 


LiFlA 
LiF2B 


9 9+1.2 
^•^-0.8 

< 2.4 (2a) 


49 
25 


GHRS 


NI 
NI 
NI 


1199.550 
1200.223 
1200.710 


2.192 
2.015 
1.713 


G160M 
G160M 
G160M 


27.6lt^ 
20.2tll 

ii.ottl 


17 
17 
16 


Mg II 
Mg II 


2796.352 
2803.531 


3.236 
2.933 


Ech-A 
Ech-A 


551^ 
37.91- 


12 
13 


Si II 


1193.290 


2.775 


G160M 


24.711;° 


16 


Fell 


2382.765 


2.882 


Ech-A 


20.512J 


12 
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Table 4 — Continued 



Species logA/^ Channel Wx [mAY S/N^ 



'^Central wavelength of each transition from Mor- 
ton 2001. 

Adopted product of wavelength and oscillator 
strength for each transition, taken from Morton 
2001. 

■^Measured equivalent width and la uncertainties 
(in mA). The uncertainties contain contributions 
from both statistical and systematic error sources. 

''Empirically estimated signal-to-noise ratios (per 
unbinned pixel) for continuum regions near the ab- 
sorption lines of interest. 

^The N III absorption contains a substantial con- 
tribution from Si II at 989.873 (fs -^21 km s"^ from 
the center of the N III line) . The equivalent width 

quoted may be dominated by the latter, but rep- 
resents a firm upper limit for the amount of N III 
along the HZ 43A sight line. 

^This transition is a blend of three O I lines, the 
sum of whose oscillator strengths are given here. 
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Table 5. Curve of Growth Column Densities 



Species 


log N ± 2a'' 


6±2cT^ 


Alo 


o 


CII 




d e 


+1. 


,20 


cm 


19 c;q+o.09 


d 


+0. 


,15 


N I 


iO.Oi_g 05 


3.01^:^ 






Nil 


13.621°;}? 


d 


+0. 


,32 


NIII 


< 12.98l°J^ 


d 


+0. 


,07 


01 




9 79+0.33 
^- '^-0.26 






VI 


12 qi +0.08 


oof 






Si II 


1 9 Qc;+0.18 




+0. 


,21 


Ar I 






+0. 


,00 


Mg II 


12.41±0.02 


2.8±0.3 






Fell 


12.17±0.04 


2.78±1.15 







'^Log of column density derived from a curve 
of growth fit, with estimated 2<t uncertainties. 
For most species not enough transitions were 
observed to independently fit a curve of growth 
and determine a 6- value. In these cases we have 
adopted the 6-value of 2.72toil (2a) from O I. 

'^Doppler parameter in km s~^ with esti- 
mated 2a uncertainties. 

'^For those species for which we have adopted 
the 6-value from O I we give the amount by 
which the curve-of-growth-corrected column 
density given in column 2 of this table exceeds 
the average integrated apparent column den- 
sity. 

'^We have adopted the 5- value of 2.721^;^^ 
{2a) from O I to estimate the total column 
density of these species. 

'^For C II the b-value is probably larger, and 
logA^ correspondingly smaller. Column 4 gives 
the maximum correction for this effect. 

^For O VI we assume that the instrument 
fuUy-resolves the absorption profile and adopt 
the weighted average of the apparent column 
density integrations from the two LiF channels. 
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Table 6. HZ 43 Column Densities and Ratios'^ 





Hi 


Di 


Oi 


Nl 


Log N 


17.93 ± 0.06 


1 q 1 r +0.040 


14.49 ± 0.08 


13.51 ± 0.06 


Ratios 


Dl/Hl 
1.66 ±0.28 X 10-^ 


Oi/Hi 
3.63 ± 0.84 X 10-^ 


Nl/Hl 
3.80 ± 0.74 X 10-^ 


Dl/Ol 
4.57 ±0.96 X 10-2 


Dl/Nl 
4.37 ± 0.74 X 10-^ 



^all uncertainties are 2 a 



